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^ ■ Abstract 

CN ' We consider a neutrino dark energy model and study its implications for the neutrino su- 
perluminality reported recently by the OPERA collaboration. In our model the derivative 



couplings of the neutrino to the dark energy scalar result in a Lorentz violation in the neu- 
^ [ trino sector. Furthermore, the coupling of the dark energy scalar field to the stress tensor of 
^ , the Earth automatically leads to a nontrivial radial profile for the scalar which in turn yields 
a terrestrial neutrino v — c far above its value in interstellar space, so as to be simultaneously 
compatible OPERA, MINOS and SN1987A data upon fitting a single parameter. 
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1 Introduction 



Recently the OPERA experiment announced the arrival of 16,111 Swiss neutrinos, 61 
nanoseconds ahead of schedule PJ. This complimented a similar, but less confident, an- 
nouncement concerning less energetic American neutrinos four years ago pj. Although the 
former of these announcements is claimed to be a 6 sigma event, as reviewed in Ref. [3] these 
claims have been met with a degree of skepticism because they need to be reconciled with a 
series of yet more solid observations: 

1) Neutrinos from SN1987A [11[5], after a 160,000 year trip, arrived on Earth at most 3 hours 
before the corresponding light [6] . The arrival time was essentially independent of the energy, 
with all neutrinos arriving within a 12.4 second window. 

2) The non-trivial energy dependence of the survival probability of electron reactor neutrinos 
at KamLAND requires that the relevant Lorentz- violating couplings must be neutrino flavor- 
independent to within a precision of about 10~^° [7]. 

Naively the first of these observations is the most inconsistent with OPERA and MINOS's 
claims. OPERA found a fractional superluminality (v-c) /c of 2.5±0.28±0.3 x 10-^ MINOS 
found 5.1 ± 2.9 x 10^^ while SN1987A establishes a maximal fractional superluminality of 
2 X 10"^. One may seek consolation in the fact that SN1987A neutrinos are electron neutrinos 
whereas the others are muon neutrinos, but the second observation above makes this point 
irrelevant. Another difference is that the SN1987A neutrinos are in fact mostly antineutrinos, 
but no model has yet come forward which successfully exploits this difference. 

A perhaps more useful distinction is the energy of these neutrinos. The SN1987A neu- 
trinos lie between 7 and 39 MeV, while MINOS neutrinos peak at about 3 GeV with a 
long tail and OPERA neutrinos are largely between 10 and 50 GeV. Given the large errors, 
one may wish to simply throw out the MINOS results and fit the fractional superluminality 
with a powerlaw E°' where a > 1.3. A steep enough power law {a > 2) can be consistent 
with the nearly simultaneous arrival time of SN1987A neutrinos, as the superluminality may 
simply not have been observed. However OPERA also analyzed the energy-dependence of 
its neutrinos, separating them into two groups with average energies differing by a factor 
of 3 and found a = 0.1 ± 0.6, which is difficult but not impossible to reconcile with the 
aforementioned consistency condition. The next year of OPERA data should allow a more 
definitive conclusion on the plausibility of this scenario. 
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Even if next year's data confirms OPERA's superluminality claim, the possibility of a 
systematic error at OPERA remains. We will not attempt to determine whether such an 
error exists, but merely note that if future T2K and MINOS runs confirm OPERA's result 
then it must be taken seriously. With this motivation we will suspend our disbelief and 
describe a variation of the neutrino dark energy model of Ref. [9] which renders it consistent 
with neutrino data from SN1987A. 

2 Neutrino Coupling and Dispersion Relations 

To change the velocity of a neutrino, one needs to change its dispersion relation. This 
change necessarily violates Lorentz symmetry. We will consider spontaneous violations of 
Lorentz symmetry, which arise by adding terms to the Lagrangian which couple a new field 
to neutrino bilinears. The couplings and the new field are similar to those introduced in 
the neutrino dark energy model of Ref. and to a large extent to those in the earlier 
models of Refs. [TU] , although we do not demand that the additional field actually provides 
the observed dark energy. The field acquires a VEV due to interactions with the Earth, 
which spontaneously breaks the Lorentz symmetry. In Refs. [131 Ull US] such models were 
constructed in which the new fields introduced were respectively a symmetric tensor, a 
vector and a scalar. We will consider a scalar field 11, which as explained below will have 
the advantage that it requires the tuning of only a single parameter. The fact that our 
couplings resemble those which arise in the neutrino dark energy scenario of Ref. [9] yields 
a cosmological justification for the exclusiveness of these terms to the neutrino sector. 

In an effective field theory setting, it is sufficient to consider the operators of lowest 
dimension which preserve Lorentz-invariance. The terms with no derivatives can be absorbed 
into redefinitions of the fields and parameters of the effective theory. Strong upper bounds 
on these terms arise, for example, on masses from beta decay and as a result these terms 
will be negligible at OPERA energies. We will be interested only in terms which cannot be 
reabsorbed into other terms via field definitions up to terms without derivatives. In general 
the modifications of the dispersion relation can be linear or quadratic in the new couplings, 
we will eventually restrict our restriction to the linear modifications and so to the coupling 
terms which lead to linear modifications. With all of these criteria, we are left with 

A£ = i (za^z/9V + ic^.u^d'^iy - d^.^i^^^^d^dPu) (2.1) 

where a, c and d are tensors constructed from derivatives of 11. In our dimensional analysis 
scheme, in which the coefficients are constructed from a scalar field of dimension [m], the d 
in the last term is of higher dimensionality and so need not be considered, but we will keep it 



2 



during this section for illustration as a linear energy dependence of neutrino superluminality 
is excluded at OPERA by less than two sigmas. 

Clearly these derivative terms are nontrivial only if 11 is not constant. OPERA and 
MINOS demand superluminality on Earth, while SN1987A neutrinos demand much less 
superluminality in space, therefore these derivatives need to be localized on Earth. The 
simplest possibility would be if 11 were localized on Earth and vanishes in space, while dkH 
vanishes everywhere. However it is easy to see that if 11 begins with a constant value, the 
spatial gradients between here and SN1987A will dominate over the temporal gradient in 
less than a mere 160,000 years, much less than the age of the Earth. Therefore such field 
configurations are not logically consistent. 

The next simplest possibility is that 11 is negligible but the gradient of 11 lies along the 
Earth's radial direction r. In Sec. E] we will describe a prototypical example of a scalar 
field with the desired behavior. For now, it will simply be relevant that the only nonzero 
derivatives, and so the only nonzero components of the tensors a, c and d, will be those with 
only r indices or an even number of identical tangential indices. Then the dispersion relation 
is 

E = + M'"^ (2.2) 

with ith spatial component 

P- = Pi + CijPj + dijkPjPk M' = m + ajpj (2.3) 

where sums over repeated indices are understood. 

Notice that as in the model of Ref. [H] and unlike that of Ref. [13], the dispersion relation 
is anisotropic. In particular, the velocity of a neutrino at the point p traveling at an angle 
6 with respect to the radial direction is 

2 

~ — cos^ 6 + (c^^ cos^ 9 + Cyy sin^ 9) + 2E cos 9 [d^^^ cos^ 9 + ?>dxyy sin^ 9) (2.4) 

where, without loss of generahty, we have chosen our coordinates such that, at p, ?/ = 2; = 
while the neutrino's velocity is in the x — y plane. 

As mentioned above in the sequel we will ignore the last term as it is created by a higher 
dimensional operator in our effective field theory description. Notice that Cxx and ax only 
occur in the combination Cxx + o^/2 at this precision. The terms which we have omitted are 
strongly suppressed, by factors of the superluminality fraction or even by the ratio of the 
neutrino energy to its rest mass, and so the identification of further terms or even the separate 
identification of Cxx and will need to await a much more precise or qualitatively different 
experimental setting. As we can then not hope to experimentally distinguish between the 
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effects of Cxx and a^., we will simply neglect in what follows. Working with the effective 
lagrangian with a SUl{2) x t/y(l) symmetry, the second operator in Eq.(2.1) is given by 



AC = -b{d''d''n){H^L)-f^d^{LH) 



(2.5) 



were we have defined the constant b 



(2.6) 



c, 
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and the H is the Higgs doublet of the standard model. In Ref. [T^ the author considers 
a term which couples the neutrino kinetic tensor to (9^11(9,^11. However this term is of one 
energy dimension greater than (12. 5p . and so is suppressed according to the usual logic of 
effective field theories. 

At this point one may already fit to the OPERA data to find c on Earth. Instead we 
will find it convenient to first describe our simple scalar model which determines c in terms 
of the derivatives of a scalar field in Sec. [31 We will then fit this parameter to the OPERA 
result in Sec. HI 



A crucial difference between neutrinos observed by OPERA and those emitted by SN1987A 
is that the first traveled within the Earth, while most of the journey of the later was in in- 
terstellar space. Thus the discrepancy between these observations can be accounted for if 
the velocity of a neutrino is higher within the Earth than in space. This can be arranged 
in a number of inequivalent ways, by considering models in which the neutrino propagator 
is modified by a kinetic coupling of the neutrino to a field which obtains a VEV, if this 
messenger field is coupled to the Earth. In Ref. [IS] , for example, the authors proposed that 
this spatial dependence can emerge in type IIB string theory model. 

No such coincidence is required in models in which the messenger is coupled directly to 
baryon density or to the background stress tensor in such a way that it acquires a classical 
expectation value concentrated near massive objects. So long as this classical field drops off 
sufficiently quickly from the sources, it will affect terrestrial neutrinos and not appreciably 
affect supernova neutrinos. However it is important that the field drops off sufficiently 
quickly so that the Earth's effects dominate preferably over those of the Sun and certainly 
over those of the center of the galaxy. In Ref. [13] this was achieved by adding a spin two 
field whose inverse mass is fixed by hand to be roughly the inverse radius of the Earth 
(and necessarily less than the inverse distance to the Sun), while the coupling was chosen 
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to yield the OPERA superluminality. In this note we will present an alternative model in 
which the coupling is still chosen to agree with the OPERA result, but the effects of the Sun 
and Galaxy are suppressed not by tuning another parameter, but simply by the derivative 
structure of our coupling to the neutrinos. 

We consider a model with the neutrino dark energy term f l2.5p coupling the neutrino to 
a scalar field with or without a minimal Galileon coupling y/fj given by the boundary DGP 
model HE] 

£n = -^s^ns^n - |(nn)(9,n)2 + a^/s^ut. (3.1) 

The coefficient a, which can be taken to be zero, is a parameter of dimension [l^] which 
parameterizes the nonlinear Galileon interaction. T is the trace of the stress tensor of all of 
the matter, except for the scalar. The coupling of 11 to the stress tensor could in principle 
lead to fifth forces beyond experimental bounds, as described in a very similar setting in 
Ref. [in], however a quick calculation shows that only the product of the coefficient of this 
coupling and the coefficient b in (12. 5 p appears in the neutrino v — c, therefore any reduction of 
the coefficient of the stress tensor coupling which may be mandated by fifth force constraints 
can be compensated by an opposite rescaling of b. While the Lagrangian itself has higher 
derivative terms, terms in the equations of motion have at most two derivatives acting on 
each n, which allows the existence of ghost-free solutions such as that which we will use. 

The Galileon interaction term is useful because it reduces short distance singularities, 
via the Vainshtein mechanism [20], at least in the presence of spherically symmetric stress 
tensor sources. More precisely, for an external source of mass M there will be a distance 
scale 



i?=^^^%l^(4aM)^/^ (3.2) 
V3 

at which the behavior of the 11 field changes. We can choose the Galileon coupling a such that 
this distance is either larger than or smaller than the radius of the Earth r. For concreteness, 
in the rest of this note, we will chose a to be sufficiently small so that R << r, and so we 
will set a = in the rest of this section. A small value of a is useful for, among other things, 
avoiding the potential formation of closed timelike curves [21]. However a will again become 
important for distance scales smaller than R, for example it may be invoked for neutrino 
phenomenology in neutron stars and, depending on its value, in the cores of massive stars. 

Let us now calculate a static, spherically symmetric field configuration 11 (r) in the pres- 
ence of a nonrelativistic matter source with density p(r). The 11 equation of motion is 
then 

- A^/SttGnP = V^n = ^dr{r^drU). (3.3) 
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This is easily integrated to yield 



/•r=r rXr^ 

r^drli = -A^^-nGN / f'^pdf = -\ — -M{r) (3.4) 
Jf=o V vr 

where M(r) is the density contained in the object up to a radius r, and we have fixed the 
constant of integration by imposing that 11 be differentiable at the origin. 

One can now easily find the second derivatives of 11, which appear in the coupling (12.51) . 
Choosing our coordinates such that y = z = at a given point 



9jn = ajn = V = -./!^^ (3.5) 



din = v2n-9jn-52n = -4v/3^(^p(r)-^^^ 

while terms with mixed derivatives vanish. In particular, in the Earth's crust 



(3.6) 



dlYi = 4v/3^(^po - P), dlU = dlU = -4^^Po (3.7) 

where p and po are respectively the average densities of the crust and of the Earth as a 
whole. More generally, Eq. (13. 7p may be applied to a point at any radius r in a background 
with an arbitrary spherically symmetric density profile if one identifies p with the density at 
radius r and po with the average density at radii less than r. 



4 Fitting 

4.1 OPERA, MINOS and SN1987A 

The equation (12. 5p modifies the dispersion relation for the neutrinos, allowing superluminal 
propagation. In the language of Ref. |T3] this corresponds to a modification of the effective 
metric in which the neutrinos propagate. The fractional superluminal neutrino velocity at 
an angle 6 with respect to the radial direction is given by inserting Eq. (12.61) into (12. 4p 

e = ^ = -^^idlU{r) cos2(^) + 9jn(r) sm^{e)), (4.1) 

where (f) is the Higgs VEV. In the case of OPERA and MINOS, cos^(6') ~ .01 and so we 
will ignore the first term. Using Eq. (13. 7p one finds 



26W ^(t;) Vo ~ (800 eVfh (4.2) 
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where we have used po ~ 5.5 gm/cm'^. And so the OPERA result e ~ 3 x 10 ^ is reproduced 
if 

h ~ 3^ . (4.3) 

(6keV)5 ^ ^ 

Of course, this fit is only reasonable if the main source of superluminality is the Earth. 
In other words, the main contribution to 11 must arise from its coupling to the Earth, and 
not for example the Sun or the matter in our galaxy. In the case of a distance object, the 
p term in Eq. fl3.6p vanishes, while the other terms are simply the ratio of the total mass 
to the distance cubed. Therefore, at the surface of the Earth, the contribution of the Sun 
(with respect to that of the Earth) is suppressed by a factor of 10~^ and that of the galaxy 
by 10~^^, thus it is reasonable to ignore these contributions. In the case of the model in 
Ref. ^13] the corresponding M/R dependence would lead to a dominant contribution from 
the Sun and galaxy, and so these contributions are eliminated by hand by choosing the mass 
of their new field to be of order the inverse radius of the Earth. 

And so what about neutrinos from SN1987A? The fractional superluminality for these 
neutrinos needs to be less than 10~^. This is easily satisfied, as one can see by directly 
calculating the change in arrival time, or else simply noting that nearly all of the trip is in 
interstellar space, where e is of order 10~^°. In fact, examining Eq. (13. 6p more closely, one 
sees that during the first part of their trip the neutrinos are actually subluminal. 



4.2 KamLAND 

Fitting the OPERA and MINOS results is easy, one need merely add a Lorentz-violating 
term to the muon neutrino sector. The SN1987A neutrinos are then not a problem, so long 
as the electron neutrino sector is unaltered. This is also consistent with LEP synchrotron 
limits, as these concern electrons, if one ignores fiavor mixings. The limits on muon Lorentz 
violation placed by gamma rays arriving from blazars [22] are violated, but since our scalar 
is only appreciable during a small portion of this trip, the relevance of these bounds to 
OPERA physics is unclear. However as we have alluded earlier, the biggest problem arises 
from KamLAND's reactor neutrino oscillation data. This suggests that certain terms which 
violate Lorentz-invariance in the muon neutrino sector must be equal in the electron neutrino 
sector. No such analysis has been carried out, however, for every possible term in our low 
energy effective description. In this subsection we will comment on the possibility of evading 
these bounds within the validity of the effective theory. 

KamLAND lies between 55 nuclear reactors, which supply it with antineutrinos. The 
electron neutrinos oscillate into muon neutrinos. If the energy difference between an electron 
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and muon neutrino is AE, the probability that a given electron neutrino is still an electron 
neutrino after traveling a distance L is 

l-sin^(2^e,)sin2(^^) (4.4) 

where 6efj, is the mixing angle between these two flavors. Using the dispersion relation f|2.2p 
one readily flnds, for an energy E neutrino, that the latter phase is 

^ ( W + + - - ) ^ ^ + 2 ^'yy ' 'yy^ ^ ^+ 

(2 (41 - dit) cos' e + 3 (4^J, - di^l) sin' e) E' cos' 6) (4.5) 

where dm' is the difference between the squared masses of the two relevant neutrino mass 
eigenstates and the e and /i superscripts denote the Lorentz-violating neutrino couplings to 
the two flavors respectively. Now inserting coefficients 2c+a' which are sufficiently large to fit 
the OPERA data, together with an energy E = 5 MeV typical of reactor neutrinos seen above 
the background at KamLAND, one flnds that neutrinos oscillate about every 5 nanometers. 
This is much smaller than the resolution at KamLAND, and so KamLAND would in this 
case flnd an energy-independent survival probability. This is in grave contradiction with 
KamLAND's results in Ref. [7] which present two full averaged neutrino oscillations. Higher 
order terms, such as the d term, lead to an even shorter wavelength but with an E dependence 
which cannot cancel that of previous terms. Therefore a cancellation in such terms, in order 
to be consistent with KamLAND's results, must be imposed order by order in E, and in fact 
also order in order in cos(6'), which is small but still larger than the ratio of the energy to 
the neutrino mass. Such a cancellation is not possible with the terms in the lowest orders of 
the effective action considered in this note. 



5 Experimental Signatures 

Once one has included the interaction term fl2.5p to explain superluminal neutrino ve- 
locities here on Earth, one needs to deal with its consequences throughout our universe. In 
particular, it may become appreciable inside of various astrophysical objects. In the case of 
a spherically symmetric object, Eq. (13. 7p tells us that the super luminality at a radius r only 
depends on the density p at that radius and the average density at lower radii po- This may 
seem peculiar, since there is no Gauss' law for 11. Indeed, any choice of boundary conditions 
for n implies that its value in the core of an object depends on the p proflle at higher radii. 
However in our model, and not in that of Refs. [T3| [H] , the super luminality only depends 
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on the second derivatives of 11, which is independent of the constants of integration and so 
of these boundary conditions. This freedom in choosing boundary conditions is necessary 
for many cosmological apphcations of the 11 field, such as the dark energy model of Ref. [9] 
and Galileon cosmologies. 

In the core of a spherically symmetric object, in which the density is approximately 
constant, the superluminality is easily approximated. Up to geometrical factors whose mag- 
nitude and sign depend on the direction in which the neutrino travels, the fractional super- 
luminosity of neutrinos is simply given by the ratio of the OPERA result to the ratio of the 
average density of the Earth to that of the object's core. For example in the case of the core 
of our Sun, whose density varies between 4 and 30 times that of the Earth, one expects a 
neutrino fractional subluminosity or super luminosity of order between 10~^ and 10~^. This 
change in the neutrino's dispersion relation affects its propagator and so may in principle 
have observable consequences on the various fusion processes occurring in our solar core, 
perhaps allowing for an exclusion or verification of our model. This is in contrast with the 
model of Ref. [TSJ, in which neutrino superluminality is of order 1 at the surface of the Sun 
and persists at greater depths including the outer core, where a significant amount of fusion 
occurs. 

However, even in our model the superluminality will be of order 1 whenever the density 
reaches about 2 x 10^ g/cim?. This level is easily surpassed for example in the cores of older 
massive stars which fuse carbon or heavier elements Ref. [23], and is even reached in white 
dwarves. Thus in these situations one may expect drastic departures from the standard 
Lorentz-invariant quantum field theory predictions for amplitudes and decay rates involving 
the neutrino propagator, potentially affecting for example models of the helium flash. 

Another crucial distinction between our model and that of Ref. [13] is that our neutrino 
velocity is direction-dependent, radially traveling neutrinos on the surface of the Earth are 
subluminal. For shallow angles such as those of OPERA, MINOS and T2K the neutrino 
superluminality in our model is reasonably independent of the angle, there is only a correction 
of order the angle (f — squared. On the other hand, in the second model of Ref. [H], 
the superluminality is proportional to the cosine of the angle with respect to a preferred 
direction, which given the symmetries of the problem will likely correspond to the Earth's 
radial direction. As the baseline of T2K is half that of the others, this cosine is only one half 
of its value in OPERA and MINOS, and so the second model of Ref. [H] suggests that the 
neutrino (i; — c)/c at T2K is about one half its value at OPERA, corresponding to a neutrino 
arrival time only 14 nanoseconds ahead of schedule. Thus if future T2K runs indicate no 
superluminality, that model is unlikely to be ruled out. On the other hand, to distinguish 
our model from that of Ref. [T3| one needs a significantly larger angle, translating into a 
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larger baseline. Almost an order of magnitude of improvement in the timing precision of 
OPERA (at both experiments) would be necessary to distinguish the models by comparing 
superluminality at OPERA and at Fermilab's proposed 1300 km baseline experiment LBNE. 
The prospects are improved significantly in the 2000+ km proposed experiments sending 
neutrinos from CERN to the Pyhasalmi mine. 

6 Electron Superluminality? 

In addition to the experimental constraints listed in the introduction, one may add: 

3) Bounds on terrestrial electron superluminality are extremely tight. At LEP [8J, synchrotron 
radiation measurements at OPERA energies yield a maximum fractional superluminality 

{v - c)/c of 5 X 10-1^ 

Together with constaint (2), this implies that either the neutrino superluminality disappears 
before the electroweak scale, or else some conspiracy of other effects prevents it from being 
transmitted to the electron. The first scenario is very difficult to realize, especially in light 
of the fact that OPERA neutrino energies are already so close to the electroweak scale. For 
example in Ref . |3] the authors noted, as was shown in Ref . [TU] , that in the case of a coupling 
ly^dQU, a one-loop mixing generates the operator e^doe which leads to a superluminality for 
electrons which is suppressed by at most E"^ /My^ with respect to that of neutrinos, and so 
is 5 orders of magnitude above the LEP bound quoted above. More complicated breaking 
models, such as the coupling to sterile neutrinos |3] and supersymmetric completions 
[5], are able to push these effects onto higher derivative operators. However generically this 
merely increases the exponent of the E/Mw suppression, but since E/Mw is of order 0.1 it 
seems difficult to see how such models could be consistent with the LEP constraint. 

One might hope that there exist flavor-asymmetric Lorentz-violating terms which evade 
the neutrino oscillation constraints of observation (2). In this case one could simply declare 
that Lorentz-violation only occurs for muon and perhaps tau neutrinos, hoping that their 
mixing with electrons is sufficiently suppressed to avoid the LEP synchrotron bound. There 
are a number of reasons why this approach fails, but we have seen (Subsec. 14. 2p that our 
terms certainly do not allow such an evasion of the flavor neutrality constraints. 

This leaves the second possibility for generating the required superluminality hierarchy. 
OPERA data, together with observations (2) and (3), appears to imply an extreme flne 
tuning which cancels the electron couplings leading to superluminality. Such a flne tuning 
can be achieved, at a flxed energy scale, by tuning bare values for the electron analogues of the 
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coupling fl2.5p and indeed it has been claimed that a cancellation is realized automatically 
in some models [HI |T2]. Of course the electron coupling will in turn generate one-loop 
corrections to the neutrino coupling, and so a proper fit to the data in principle requires 
one to fit these two bare couplings simultaneously, and one may also wish to include such 
couplings for the other charged leptons, to satisfy for example Lorentz-violation bounds in 
the muon sector. While the (provisional) resolution of hierarchy problems via the fine-tuning 
of bare parameters is certainly familiar and perhaps even necessary in particle physics, in 
our opinion the need for this fine-tuning is among the most compelling reasons to suspect a 
systematic error in the OPERA data. 
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